1.1 Contributions

Please note that authorship is alphabetical. Contributions are listed below - see github for details and who to blame for what :-).

Ben Anderson (b.anderson@soton.ac.uk @dataknut)

1.2 Code

• https://github.com/CfSOtago/airQual/commits/master

1.3 Citation

If you wish to refer to any of the material from this report please cite as:

• Anderson, B., (2019) Air Quality in New Zealand: Exploring “official” data , University of Southampton: Southampton, UK.

Report circulation:

• Public.

Report purpose:

• to explore official New Zealand Air Quality data

This work has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 700386 (SPATIALEC).

This work is (c) 2019 the University of Southampton.

5.1 PM10

skimr::skim(pm10dt)

## Skim summary statistics
##  n obs: 209964 
##  n variables: 10 
## 
## ── Variable type:character ────────────────────────────────────────────────────────────────────────────────
##      variable missing complete      n min max empty n_unique
##       council       0   209964 209964   2   5     0       15
##  council.site       0   209964 209964   8  34     0       93
##          date       0   209964 209964  19  19     0     4382
##        method       0   209964 209964   3   4     0        2
##          site       0   209964 209964   4  28     0       93
## 
## ── Variable type:Date ─────────────────────────────────────────────────────────────────────────────────────
##  variable missing complete      n        min        max     median
##   ba_date       0   209964 209964 2006-01-02 2017-12-31 2011-12-27
##  n_unique
##      4382
## 
## ── Variable type:logical ──────────────────────────────────────────────────────────────────────────────────
##            variable missing complete      n mean
##   complete_for_mean    1152   208812 209964 0.74
##  complete_for_trend    1152   208812 209964 0.64
##       complete_year   74218   135746 209964 1   
##                              count
##  TRU: 154552, FAL: 54260, NA: 1152
##  TRU: 134660, FAL: 74152, NA: 1152
##             TRU: 135746, NA: 74218
## 
## ── Variable type:numeric ──────────────────────────────────────────────────────────────────────────────────
##  variable missing complete      n  mean    sd  p0 p25 p50   p75   p100
##      pm10       0   209964 209964 17.12 12.83 -10 9.9  14 20.22 433.75
##      hist
##  ▇▁▁▁▁▁▁▁

# looks like daily data with gaps
t <- pm10dt[, .(nObs = .N), keyby = .(council, site)]

kableExtra::kable(t, caption = "N obs at sites (PM10)") %>% kable_styling()

5.2 PM2.5

skimr::skim(pm2.5dt)

## Skim summary statistics
##  n obs: 33750 
##  n variables: 10 
## 
## ── Variable type:character ────────────────────────────────────────────────────────────────────────────────
##      variable missing complete     n min max empty n_unique
##       council       0    33750 33750   2   4     0        9
##  council.site       0    33750 33750   7  26     0       35
##          date       0    33750 33750  10  10     0     3617
##        method       0    33750 33750   3   3     0        1
##          site       0    33750 33750   4  22     0       35
## 
## ── Variable type:Date ─────────────────────────────────────────────────────────────────────────────────────
##  variable missing complete     n        min        max     median n_unique
##   ba_date       0    33750 33750 2008-01-01 2017-11-25 2014-03-05     3617
## 
## ── Variable type:logical ──────────────────────────────────────────────────────────────────────────────────
##            variable missing complete     n mean
##   complete_for_mean       0    33750 33750 0.56
##  complete_for_trend       0    33750 33750 0.32
##       complete_year   14594    19156 33750 1   
##                          count
##  TRU: 18763, FAL: 14987, NA: 0
##  FAL: 22856, TRU: 10894, NA: 0
##          TRU: 19156, NA: 14594
## 
## ── Variable type:numeric ──────────────────────────────────────────────────────────────────────────────────
##  variable missing complete     n mean   sd    p0 p25  p50  p75 p100
##     pm2_5       0    33750 33750 8.19 7.89 -0.43 4.1 5.79 8.69  180
##      hist
##  ▇▁▁▁▁▁▁▁

# looks like daily data with gaps
t <- pm2.5dt[, .(nObs = .N), keyby = .(council, site)]

kableExtra::kable(t, caption = "N obs at sites (PM2.5)") %>% kable_styling()